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RESUMEN 

Combinando modelos de smtesis de poblaciones estelares con simulaciones Montecarlo de los efectos estocasticos 
en el numero de estrellas que ocupan estadios evolutivos escasamente poblados en el DHR, se concluye que 
la dispersion observada en la magnitud y colores fotometricos de cumulos estelares en la NGM y NGC 7252 
pueden entenderse en el marco de la teorfa actual de evolution estelar. El uso de bibliotecas de espectros 
estelares de alta resolution en los modelos de smtesis de poblaciones mejora considerablemente la calidad de 
los ajustes a las distribuciones espectrales de energfa observadas, y permite asignar edades espectroscopicas a 
poblaciones estelares con mayor confiabilidad que con los modelos de baja resolution espectral. 

ABSTRACT 

Combining population synthesis models with simple Montecarlo simulations of stochastic effects in the number 
of stars occupying sparsely populated stellar evolutionary phases in the HRD, I show that the scatter observed 
in the photometric magnitudes and colors of LMC and NGC 7252 star clusters can be understood in the 
framework of current stellar evolution theory. The use of a high resolution stellar spectral atlas in population 
synthesis models improves considerably the quality of the fits to observed galaxy SEDs, making the assignment 
of a spectroscopic age to stellar populations more reliable than with low spectral resolution models. 
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1. INTRODUCTION 

Evolutionary population synthesis has become an 
important tool to study the stellar population con- 
tent and its evolution over cosmological time scales 
in star cluster and galaxies. Bruzual (2001, 2002a) 
has presented a detailed description of the popu- 
lation synthesis technique, including main results, 
comparison with observations, and sources of uncer- 
tainties. For reasons of space, in this paper I include 
a very brief summary of population synthesis and 
an update on recent developments since the previ- 
ous papers. References to complementary work are 
provided. 

2. POPULATION SYNTHESIS 

The current theory of stellar evolution makes de- 
tailed predictions concerning the time dependence of 
the bolometric luminosity and effective temperature 
of a star of a given mass and chemical composition. 
Complete sets of evolutionary tracks for stars of a 
wide mass range and various initial metal contents 
covering all significant stellar phases are available in 
the literature (e.g. Fagotto et al. 1994a, b, c; Gi- 
rardi et al. 2000). Assuming an initial mass func- 
tion (IMF, e.g. Salpeter 1955), we can compute the 



number of stars of a given mass born at time zero, 
and then follow the evolution of this population in 
the HRD using a specific set of evolutionary tracks. 
Knowledge of the spectral energy distribution (SED) 
at each position in the HRD visited by the stars 
during their evolution, allows us to compute the in- 
tegrated SED for this initial-burst or simple stellar 
population (SSP) as a function of time. By means 
of a convolution integral (Bruzual & Chariot 1993) 
the evolution of the SED can be computed for an 
arbitrary star formation rate (SFR) and a chemical 
enrichment law. 

3. STELLAR LIBRARIES 

Several stellar spectral libraries for solar metal- 
licity (Zq) stars are currently available. The Pick- 
les (1998) atlas provides good coverage of the HRD 
for stars of this metallicity at medium spectral res- 
olution. Le Borgne et al. (2002) have compiled an 
equally comprehensive atlas at 3A spectral resolu- 
tion (lA sampling) which is complete for stars of so- 
lar metallicity, but includes a large number of spec- 
tra of non-solar metallicity stars. On the theoretical 
side, the compilations by Lejeune et al. (1997, 1998) 
and Westera et al. (2002) provide libraries of model 
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Fig. 1. SED of a Z = Z Q BC02 SSP model at 
10, 100, and 1000 Myr computed with different stellar 
spectral libraries: (light gray) ~20A resolution (Lejeune 
et al. 1997), (medium gray) sampling interval of 5 A 
pixel" 1 and resolution of A/ A A « 500 (Pickles 1998), 
and (black) sampling interval of 1 A pixel -1 and resolu- 
tion of A/ A A ?a 1700 (Le Borgne et al. 2002). The SEDs 
have been shifted in the vertical direction for clarity. 



atmospheres for stars of various metallicities but at 
«20A spectral resolution. The last two libraries are 
largely based on the Kurucz (1995) series of model 
stellar atmospheres. Figs 1 and 2 show the time 
evolution of the SED for a Z Q SSP model computed 
for these three stellar libraries in increasing order 
of spectral resolution. Throughout this paper I use 
the x = 1.35 Salpeter (1955) IMF, and for the lower 
and upper mass limits of star formation, I assume 
to/ = 0.09 and m u — 125 M Q , respectively. The 
set of population synthesis models that use the Le 
Borgne et al. (2002) high resolution stellar spectral 
library will be discussed by Bruzual & Chariot (2002, 
BC02 hereafter). All the models shown in Figs 1 
and 2 have the same behavior in what respects to 
integrated photometric properties of the stellar pop- 
ulation, e.g. luminosity and color evolution. The 
highest resolution model reproduces in much greater 
detail the spectral features of integrated stellar pop- 
ulations, observed, for instance, with the high reso- 
lution spectrograph available in HST and in the new 
generation of ground based high performance large 
optical telescopes (Keck, VLT, Sloan). See for ex- 
ample Figs 1 and 2 of Bruzual (2002c) and §5 below. 
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Fig. 2. SED of a Z = Z e BC02 SSP model at 
2, 4, and 8 Gyr computed with different stellar spec- 
tral libraries: (light gray) «20A resolution (Lejeune et 
al. 1997), (medium gray) sampling interval of 5 A 
pixel" 1 and resolution of A/ A A « 500 (Pickles 1998), 
and (black) sampling interval of 1 A pixel -1 and resolu- 
tion of A/ A A « 1700 (Le Borgne et al. 2002). 

4. STAR CLUSTER COLORS 

The differences between the observed integrated 
colors of star clusters and the colors predicted for 
simple stellar populations by population synthesis 
models are, in general, larger than allowed by obser- 
vational errors. This is particularly true for interme- 
diate age clusters in optical-IR colors, like V — K, 
but to a lesser extent the scatter is also large in 
B — V . In Fig 3 I compare the observed B — V and 
V — K colors of LMC globular cluster (van den Bergh 
1981; Persson et al. 1983, respectively) and the star 
clusters in NGC 7252 (Maraston et al. 2001) with 
model predictions. The lines in the top frame repre- 
sent the evolution in this plane of BC02 SSP mod- 
els computed using the Padova evolutionary tracks 
for Z = 0.2 x Z Q , 0.4 x Z Q , and Z Q (Fagotto et 
al. 1994a,b,c), the Salpeter (1955) IMF, and the 
Lejeune et al. (1997) stellar atlas. The * symbols 
along the Z — 0.4 x Z Q line mark the model col- 
ors at the age indicated by the labels, and can be 
used to roughly date the clusters. The small dots 
in the bottom frame of Fig 3 indicate the B — V 
and V — K colors resulting from different simula- 
tions by Bruzual (2002b), who has followed the San- 
tos & Frogel (1997) Montecarlo technique to show 
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Padova (1994), Z = 0.008, Salpeter IMF 
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Fig. 3. V — K vs. B — V color diagram. The large 
symbols represent LMC globular clusters discriminated 
according to their SWB class. The points with error bars 
correspond to the star clusters in NGC 7252. The lines 
represent the evolution of BC02 SSP models for various 
metallicities. The dots in the bottom frame indicate the 
model colors when stochastic fluctuations in the number 
of stars are taken into account (Bruzual 2002b). 



Fig. 4. K (at the distance of the LMC) vs. V - K 
diagram. The symbols represent LMC globular clusters 
discriminated according to their SWB class. The line 
represents the Z = 0.4 x Zq SSP model with no fluc- 
tuations. The dots indicate the results of simulations 
including stochastic fluctuations in the number of stars 
in the same model for a total cluster mass of 1 x 10 5 (top) 
and 1 x 10 4 M (bottom). 



that these cluster colors are consistent with model 
colors, if stochastic fluctuations in the number of 
stars in sparsely populated evolutionary stages are 
properly included into the models. 100 simulations 
were run at each of 220 time steps or isochrones, 
obtained from the Z — 0.4 x Zq SSP model shown 
in the top frame of the figure. Each simulation was 
stopped when the total cluster mass (including dead 
stars) reached 1 x 10 4 Mq. The fluctuations in the 
colors become larger as the cluster mass decreases. 
The model with no fluctuations is equivalent to an 
infinite number of stars populating the IMF. The ex- 
pected fluctuations in V — K for a 1 x 10 4 M Q cluster 
amount to almost 2 magnitudes, in close agreement 
with the range of colors observed at a given age, 
and is considerably broader than for the 1 x 10 5 M Q 
case. At the intermediate ages the models are red- 
der in V — K than the analytic IMF model because 
of a larger number of AGB stars, which appear nat- 
urally as a consequence of stochastic fluctuations in 
the IMF. The region in this plane occupied by the ob- 
servations and the cluster simulations overlap quite 
well. This comparison seems to favor a not very large 
mass for these clusters, if stochastic fluctuations are 



responsible of the color variations. This conclusion 
is supported by Fig 4 in which I plot K (at the dis- 
tance of the LMC) vs. V — K. The region covered 
by the 1 x 10 4 M Q cluster simulations matches the 
region covered by the observations very well. See 
Bruzual (2002b) and references therein for details 
and Cervino et al. (2000) for another approach to 
this problem. 

5. MODELING GALAXY SEDS 

Using the BC02 models we can derive the spec- 
troscopic age that characterizes a given observed 
SED. A model consists of a spectrum which evolves 
in time from age zero to 20 Gyr in 221 time 
steps. By minimizing x 2 = ^2[log F\(observed) — 
log F\(model)] 2 , we derive for each problem SED the 
(spectroscopic) age at which the model matches more 
closely the observed spectrum. Figs 5 to 8 compare 
the observed SEDs of 4 galaxies in the SDSS early 
data release with the best fitting BC02 model SED. 
Figs 5 and 6 correspond to relatively old stellar pop- 
ulations, which are reproduced remarkably well by a 
SSP at age 4 and 8 Gyr, respectively. On the other 
hand, Figs 7 and 8 correspond to galaxies which are 
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Fig. 5. Comparison of a galaxy SED from the SDSS 
early data release (gray line) and the 4 Gyr, Z = Zq 
BC02 model SSP computed for the Salpeter (1955) IMF 
using the Le Borgne et al. (2002) stellar atlas (black 
line) . 
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Fig. 6. Comparison of a galaxy SED from the SDSS 
early data release (gray line) and the 8 Gyr, Z — Zq 
BC02 model SSP computed for the Salpeter (1955) IMF 
using the Le Borgne et al. (2002) stellar atlas (black 
line) . 



actively forming stars, as indicated by the emission 



Fig. 7. Comparison of a galaxy SED from the SDSS 
early data release (gray line) and the 5 Gyr, Z = Zq 
BC02 model with star formation rate \]>(t) oc exp(— t/r) 
with t = 5 Gyr (black line). 
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Fig. 8. Comparison of a galaxy SED from the SDSS 
early data release (gray line) and the 10 Gyr, Z — Zq 
BC02 model with star formation rate 9(t) oc exp(— t/r) 
with t — 5 Gyr (black line). 



lines seen in their SEDs. The nebular lines emitted 
by ionized gas surrounding hot stars are not included 
in the BC02 models (see Chariot & Longhetti 2001; 
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Magris et al. 2002) but the fit to the continuum and 
the absorption lines in the observed spectra is oth- 
erwise excellent. By extending this analysis to inter- 
mediate and high rcdshift galaxies, Bruzual (2002c) 
has shown how the SEDs and/or multi band photom- 
etry of galaxies from z = to z = 4 can be ordered 
in a plausible sequence, along which the spectra of 
galaxies may evolve when passive evolution domi- 
nates the evolution of simple stellar populations. 

6. DISCUSSION 

In this paper I have briefly summarized the re- 
sults of simple numerical simulations which show 
clearly how the range of colors observed in interme- 
diate age star clusters can be understood on the basis 
of current stellar evolution theory, if we take into ac- 
count properly the stochastic fluctuations expected 
in the number of stars occupying sparsely populated 
evolutionary stages in the HRD. There does not seem 
to be any need to introduce ad-hoc assumptions in 
population synthesis models, which translate into a 
departure from our current understanding of stellar 
evolution theory, in order to explain the observed 
range of values of cluster colors and magnitudes, as 
suggested by Maraston et al. (2001). The conclusion 
that cluster masses in the range around 1 x 10 4 M Q 
are preferred, as well as the dependence of the am- 
plitude of the fluctuations in the photometric mag- 
nitude and colors on this value of the cluster mass, 
should be explored in detail. 

The use of a high resolution stellar spectral atlas 
in population synthesis models improves consider- 
ably the quality of the fits to observed galaxy SEDs, 
making the assignment of a spectroscopic age to stel- 
lar populations more reliable than with low spectral 
resolution models. The error in the age determina- 
tion depends on the quality of the observational data 
and on the uncertainties in the population synthesis 
models. In general, the relative error in the spec- 
troscopic age increases with decreasing galaxy age. 
As shown by Bruzual (2002c), it is possible to es- 
tablish a plausible sequence along which the spectra 
of galaxies may evolve when passive evolution domi- 
nates the evolution of simple stellar populations, at a 
rate that is consistent with that expected in the most 
accepted cosmological model describing our universe, 
H a = 70 km s" 1 Mpc" 1 , Q M = 0.28, Q A = 0.72, 
tu{0) = 13.75 Gyr. 



Work in progress suggests that the use of high reso- 
lution libraries for various mctallicitics will allow us 
to break the age-metallicity degeneracy and to de- 
rive a reliable star formation history from a galaxy 
SED (Mateu et al. 2002). 
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